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Abstract 



We present a theory of superconductivity in presence of a general magnetic 
structure in a form suitable for the description of complex magnetic phases 
^h| encountered in borocarbides. The theory, complemented with some details of 

c/j ' the band structure and with the magnetic phase diagram, may explain the 

nearly reentrant behaviour and the anisotropy of the upper critical field of 
HoNi2B2C. The onset of the helical magnetic order depresses superconductiv- 
ity via the reduction of the interaction between phonons and electrons caused 
by the formation of magnetic Bloch states. At mean field level, no additional 



O ' suppression of superconductivity is introduced by the incommensurability of 



the helical phase. 

PACS numbers: 74.70.Dd, 74.25.Dw, 74.25.Ha 
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0> I In 1957 V. L. Ginzburg first showed that superconductivity and long range ferromagnetic 

^ ' order compete with each other making their mutual coexistence nearly impossible [|l|. On 

O ■ the other hand in 1963 W. Baltensperger and S. Strassler pointed out that in the case of 

antiferromagnetic order the two order parameters may actually coexist 0. Two families 

Q> ! of magnetically ordered superconductors, the rare earth ternary compounds i?Rh4B4 and 

i^MogSg, were almost simultaneously discovered in the late seventies and both theoretical 
predictions proved to be correct, see e.g. Ref. for an extensive review. 

More recently the discovery of quaternary compounds of the family i?Ni2B2C, with 

^ . R = Lu, Y or rare earth element, 0J^ renewed the interest in the issue because more 

O '• complex magnetic structures were observed to coexist or to compete with superconductivity 

(for a review see Ref. [^]). For R = Er and Tm the magnetic structures coexisting with 
superconductivity are incommensurate transversely polarised spin-density waves: Tc = 11 K 

X ■ and Tn = 6.8 K for Er and T^ = 11 K and Tn = 1.5 K for Tm. DyNi2B2C, with T^ = 6 K, 

j^ ■ orders in a commensurate antiferromagnetic state at T/v = 10.6 K and it is the only boro- 

carbide compound with Tjy > T^. The case of HoNi2B2C (T^ = 8-9 K) is more complex: 
the transitions into two incommensurate magnetically ordered states, at Tj^ ~ T^fj ~ 6 K 
with wave vectors Qc = 0.91c* and Q^ = 0.55a* respectively, coincide with a deep depres- 
sion of the superconducting upper critical field 0] , while below the transition temperature 
Tjv ~ 5 K into a commensurate antiferromagnetic state, with Qaf = c*, Hc2 rapidly recovers. 
The anisotropic upper-critical-field phase diagram with the magnetic field along the sym- 
metry directions of the crystal is given in Ref. ^. Experimental data on pseudo-quaternary 

1 



borocarbides @JT31 add insight on the dramatic effects that the presence of magnetic order 
produces on superconductivity. On the other hand no difference is found in the magnetic 
properties of related superconducting and non-superconducting compounds. 

Borocarbides have a body-centered-tetragonal lattice structure with I^/mmm space 
group symmetry. In spite of their layered structure they possess three dimensional conduc- 
tion bands of mainly Ni So? character [|TI] , p]2| . The strongly anisotropic magnetic properties 
are associated with the localized 4/ electrons of the rare earths. The two types of electrons 
interact weakly via the local spin exchange on the rare earth sites. Relying on experimental 
evidence we consider the 4/-electrons system to be independent from the state (normal or 
superconducting) of the conduction electrons. Theoretically this is justified by the sepa- 
ration in the energy scales associated with the two ordering phenomena: Emo ~ ^bTm 
for magnetic order and Esc ~ (MbTc)'^/Ef for the superconductivity. Careful treatment is 
needed for magnetic structures with small q (e.g. |q| < 1/^, with ^ the coherence length) 
or originating from nesting features of the Fermi surface. In fact the opening of the small 
superconducting gap affects significantly the RKKY-type magnetic interaction in the regions 
of q-space close to zero or to a nesting vector. However, away from these special q-points, 
the structure of the magnetic interaction is related to electron-hole excitations of all ener- 
gies and is independent from the presence of superconductivity. Therefore, it is possible 
to take the magnetic properties of borocarbides from experiments or from an independent 
microscopic magnetic model and concentrate on the influence which the molecular field of 
the ordered moments has on the conducting electrons. Similar approaches have been used 
in the cases of the antiferromagnetic [pr3| , p!4| and small q helical order [^ coexisting with 



superconductivity in i?Rh4B4 and i^MoeSg. More recently the same technique has been used 



for a qualitative discussion of the properties of HoNi2B2C [p!q,p!7| . 

The aim of this work is to present a theory of superconductivity in an arbitrarily modu- 
lated exchange field capable of describing the magnetic structures encountered in the boro- 
carbides. Applied to the upper-critical-field phase diagram of HoNi2B2C, the theory re- 
produces naturally its several anomalous features and gives an explanation for its nearly 
reentrant behaviour. 

We introduce the following BCS model Hamiltonian for the conduction electrons in the 
presence of a periodic molecular field: 

H = 'Hb + 'Hh-if + Tih-b (1) 

^6 = XI ^k clck (2) 

k 

^^4/ = EAq)(Sq)-s,, (3) 

q 

^^^ = \ E y'S'^ 4ftfk,c,,. (4) 

k'^k'^k2ki 

The term Ti^ is the Hamiltonian for the free conduction band electrons, c^ (= c^^) and c^ 
(= Cko-) are respectively the creation and annihilation operators for a state with quantum 
number k = (k, a). The term Tib-Af is the exchange interaction between the spin of the local 
4/ moments (Sr,-) (with Fourier transform (Sq)) and the spin of the conduction electrons, 
Sq = J2kaa' (^k+qa' '^ o^' o^ '^kcr- At the momcut wc do uot spccify any magnetic configuration 



and we only require the periodicity of the known function (Sr-). TCb-b is the intra-band 
interaction term whose attractive part leads to superconductivity. Note that Vi^,\,^ = for 
all matrix elements not conserving spin and crystal momentum (in the form k2 + k'^ + G = 
k2 + ki, with G a reciprocal lattice vector). 

The first two terms of the total Hamiltonian are bilinear and therefore may be easily 
diagonalized for the magnetic structures encountered in the borocarbides (i.e., helices and 
SDWs). The magnetic Bloch-states obtained this way, with creation and annihilation oper- 
ators c^ (= cj^^) and c^ (= Ckr), are labelled by the momentum k and the quantum number 
r = + or — . As a convention we assume the momentum k to belong to the non-magnetic 
Brillouin zone, in order not to introduce an additional magnetic band index. Correspond- 
ingly the law of crystal momentum conservation is satisfied modulo a vector in the magnetic 
reciprocal lattice G. The magnetic reciprocal lattice may be constructed by adding to every 
non- magnetic vector G a finite set of vectors {G^}, and momentum conservation requires: 
k2-|-k'^-|-G + G^ = k2 + ki. In the new basis the Hamiltonian (|I|) reduces to7i = Tib + 'Hb-b 
with Tib and Hb-b given by eqs. (0) and (H) with all the symbols written with tildes. When 
spin-degeneracy is lifted the energy ek becomes e^. Now the Hamiltonian 7i is formally very 
similar to the usual BCS Hamiltonian. The main differences are the modified law of the 
momentum conservation and the additional k-dependence of the magnetic energy bands e^ 
and the electron-electron interaction Vj^i\i^- The mean field approximation may be applied 

to Ti. in the same way as in the non-magnetic case, via the introduction of the gap functions 
of the new magnetic eigenstates Aqm (k) corresponding to the anomalous Green functions 
(c-k+G^T'Ckr)- The function A'^ljik) is actually a matrix in the r indices in order to include 
odd and even parities. 

Some general qualitative properties of the magnetic Bloch-states and the implications 
of their use for superconductivity have been discussed in Ref. |jT^. However, numerical 



complications due to the non-trivial k-dependence introduced so far prevent the solution of 
the self-consistent gap equations in the general case. In order to proceed further the explicit 
form of the underlying magnetic structure of a particular material is needed. 

The obvious choice for the first application of this formalism is the analysis of the much 
debated issue of the almost reentrant upper critical field in HoNi2B2C. In what follows 
we concentrate on the magnetic ordered states with periodicity along the c-axis: the high 
temperature incommensurate helix (Qc = 0.91c*) and the low temperature commensurate 
antiferromagnetic state {Qaf = c*). The qualitative features of the lock-in transition are 
reproduced by theoretical models including the RKKY interaction and the crystalline electric 
field (CEF). Two such models |T^J2^ were first developed to account for the complex meta- 



magnetic phase diagram of HoNi2B2C at T = 2K |^. In particular the model in Ref. 



|T^, which includes the actual CEF states of Ho, is capable to produce a temperature 



dependent phase diagram, which is used here as input. All the parameters of this model 



have been fitted to the low temperature magnetic properties of the normal state [|19| and 
are not considered adjustable quantities in what follows. In fig. |I] we show the magnetic 
phase diagram of the model as a function of the temperature and of the magnetic field along 
the easy axis of the Ho moments, given by the crystallographic (110) direction. Input data 
for the coexistence analysis are the following calculated quantities: the underlying magnetic 
structure, the temperature dependent magnetic order parameter S{T), the magnetic energies 



and the magnetic susceptibility. 

Given the hehcal magnetic structure Sr,- = S'[acos (Q ■ Rj 
action Hamihonian in eq. (H) has the following explicit form: 



Ssin (Q ■ R,;)], the inter- 



K 



b-4f 



IS 



-Qo-'^o-'cr Ckff + Ck_|_Qo-/0"^/o- Cko- 



kcrir' 



(5) 



with cr^ = cr^ ± icTy. The corresponding magnetic states are found via simple Bogoliubov 
transformation of the type c^_^_ = u^c^-^ + fkCk+qj. where the state (k, |) mixes only with 
one other state, namely (k+Q,|), independently of the value of Q |2^. The Bogoliubov 
coefficients u^ and vy_, the new energies e^ and the new scattering matrix can be derived 
analytically in close analogy with the antiferromagnetic case. In particular we have: 
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(6) 



and: 



ek± 



Ck + CkiQ Ck — CkiQ 



A 



1 + 



4/252 



[Ck — ekiQj 



(7) 



The expression for Vj^i\i^ (not shown) is analytic as well. However, even assuming the simple 
BCS interaction potential in the non-magnetic state, it is rather complicated. It is important 
to note that the magnetic energy bands ik possess magnetic gaps only for the two pairs of 
magnetic-Bragg planes orthogonal to the c-axis at distances Q/2 and (c* - Q/2) from the F- 
point and that all the magnetic quantities differ from the non-magnetic ones only in narrow 
regions around these planes p^. The presence of only four active magnetic-Bragg planes 



for any commensurate value of the ordering vector Q make the following arguments hold for 
structures with arbitrary periodicity. If we assume that the onset of magnetic order affects 
the superconducting state as a perturbation, the only sizable components of the gap function 
matrix are Aq~ = — Aq"*" = A |]1^. After some algebraic manipulations the gap equation 
simplifies considerably and the gap function may be written as A(k, T) = (u^ — v^)A{T) 
an. This leads to: 



A(T) 



de V 



dS' 



[ut. 



- vt. 



MFS 



(27r)3 iVk'Ck' 



A (T)^(T) 

+ A2(T) 



(8) 



where JF(T) = (1 — 2nk') takes into account the occupation of the electronic states and 
we have approximated A(k, T) with A(T) in the square root. Equation (||) corresponds to 
the usual BCS self-consistent equation with an effective interaction parameter Ae(T) defined 
as the term into brackets. Ae(T) depends on the underlying magnetic state through the 
Bogoliubov coefficients and through the shape of the magnetic Fermi surface (MFS). Since 
all the anomalous magnetic k-dependencies come from the regions where the Fermi surface 
intersects the four Bragg planes, the difference AA(T) = A — Ae(T) between the actual 
electron-phonon interaction parameter (A) and the effective one may be expanded in terms 



of — —. This simplifies considerably the analysis because only very limited knowledge of 
the actual band structure is needed in order to estimate the variation AA(T) of the effective 
interaction. Therefore the high temperature interaction parameter A can be considered 
a phenomenological parameter to be determined by the transition temperature Tc. An 
inspection of the LDA band structure of HoNi2B2C |^ suggests approximate rotational 



symmetry at the intersection between the Fermi surface and the Bragg planes corresponding 
to the observed magnetic phases. With the assumption of rotational symmetry the only 
parameters we need for the band structure are the radial and vertical component of the 
Fermi velocity Vr, Vz at the intersections. To first order in the parameter — —^ for the two 
pairs of Bragg planes we have: 

AA(T) = --^^/5(T) (9) 

where kr is the radius of the intersection. We note that the two components of the Fermi 
velocity enter eq. @ independently and the perturbation expansion actually breaks down if 
Vr or Vz are much smaller than vp. In particular in the case of nesting we have Vr <^ vp and 
AA(T) is not a linear function of S(T). However the magnetic structures along the c-axis 
are not linked to nesting and may safely be treated within our perturbation expansion. We 
assume the relevant phonons to have an average energy ujd = 40meV and with the value 
Tc = 8.5K, the BCS formula gives us a value of A = 0.25. Taking into account the closeness 
of the two magnetic ordering vectors Q^^ = c* and Qc = 0.91c* we assume Vr, Vz and 
kr not to change with the lock-in. This means that we treat the incommensurate and the 
commensurate phases exactly on the same footing. 

The main features of the experimental anisotropic phase diagram of HoNi2B2C are re- 
ported in Ref. |^ which we will use in order to compare our results. The upper critical field 
may be calculated from the equation: 

h!:1''\t) = Bg^^{\{T),T) - M{T) (10) 

where Bf^^{\{T)^T) is the critical field value in the non-magnetic BCS case and M{T) 
is the magnetization in the normal state. Since the magnetic response of HoNi2B2C is 
very weak along the c-axis we may neglect the magnetization term and the upper critical 
field curve that we obtain is the upper line ((001)) in fig. |^ where we have used the value 
for AA(Oi^)/A = 0.12. The depression of the critical field is then due to the onset of the 
magnetic order altogether and not to its incommensurate nature. The depression results 
from a small but rapid decrease of the effective interaction parameter Ae(r) related to any 
helical structure, regardless to its periodicity. Furthermore the almost reentrant behaviour 
shown in fig. |^ is produced by a reduction in the value of A by only 12%. The small jump at 
the lock-in transition is due to the discontinuity of the single-site magnetization going from 
5-cell helix to the antiferromagnet shown in the lower part of fig. §. 

In order to reproduce the other features of the anisotropic upper critical field phase 
diagram is important to take into account the strong anisotropic magnetic response of the 
system and the presence of metamagnetic transitions in the same range of fields and temper- 
atures. In fig. ^is shown the curve for the external field along the (110) direction calculated 
with the M{T) response of the magnetic system. With respect to the (001) curve, the peak 



around 6K is strongly reduced by the magnetization, while for temperatures below 4K the 
two curves approach each other due to the saturation of the ordered microscopic magnetic 
moments and the subsequently reduced magnetic response of the system. The low temper- 
ature plateau is due to the transition to the metamagnetic phase AF3 (compare with fig. |l]) 
which has a large ferromagnetic component and suppresses superconductivity much more 
strongly than the low field antiferromagnet. The (100) critical field curve (not shown) has 
almost the same shape as the (110) curve for temperatures higher than 3K, but the plateau 
is reached at a slightly smaller temperature and for a slightly larger value of the magnetic 
field. This corresponds to the upper shift of the metamagnetic transition to AF3. All these 
features are in quantitative agreement with the experimental data [||. 

In addition our model explains in a natural way the fact that HoNi2B2C samples with Tc 



reduced via different techniques, i.e. with Co doping [|I0|, actually reenter the normal state 
in a temperature region around the lock-in transition T/y. The dashed line in fig. |^ is the 
upper critical field along the c-axis obtained leaving all the parameters except Tc unchanged. 

In conclusion, we derived a theory of superconductivity in a magnetically ordered back- 
ground and we applied it to the case of HoNi2B2C. We interpreted the main anomaly of its 
upper critical field via the reduction of the interaction between phonons and electrons in 
the Bloch-states of the magnetic structure. In this respect, the effect of a helical magnetic 
background on superconductivity is identical to the effect of antiferromagnetism. The helical 
case can be treated analytically independently of the periodicity and the incommensurate 
limit does not introduce additional suppression. Finally, the anisotropy of the magnetic field 
and the temperature phase diagram are well reproduced by taking into account the magnetic 
response of the material. 

We would like to thank H. Rosner for the data on the LDA band structure of HoNi2B2C. 
This work was performed under DFG Sonderforschungsbereich 463. 



FIGURES 




FIG. 1. Magnetic phase diagram for HoNi2B2C in the H-T plane obtained from the model 



discussed in Ref. |19|. The magnetic field lies along the (110) easy direction (/'). The phases have 
ferromagnetic alignment in the a6-plane and the stacking sequences along the c-axis shown in the 
figure. Inset: the RKKY interaction function v7(q). 




FIG. 2. Lower panel: size of the magnetic order parameter (S) vs. temperature. For tempera- 
tures above (below) the first order lock-in transition temperature Tjv the 115 (AF2) order parameter 
is plotted. Upper panel: Upper critical field curves Hc2(T) along the (001) and along the (110) 
directions. The dashed line the (001) curve with Tc reduced to 6.8 K 
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